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weight 4, and are functions of the invariant mass of the leptonic decay products of the 
VF-boson. 
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1. Introduction 

The measurements of inclusive semileptonic B meson decays, such as B ^ I v and 
B Xcli>, aUow a precise determination of the CKM matrix elements \Vub\ and \Vcb\- 
The latter are relevant for the study of flavor and CP violation in the quark sector (for a 
recent review see [1]). 

Total decay rates of the B meson are described by a local Operator Product Expansion 
(OPE) in inverse powers of the 6-quark mass nrn,. To leading order in 1/mf,, the total B 
meson decay rate is equivalent to the decay rate of an on-shell h quark, which can be 
calculated in perturbation theory [2]. Many authors contributed to the calculation of the 
radiative corrections to the total decay rate of 6 ulv and h ^ clD, at 0{as) [3,4] 
and 0{a^g) [5-12]. However, experimental collaborations need to impose cuts (also severe) 
on the kinematic variables. For instance, in charmless semileptonic decays, the need to 
suppress the charm background (which is ~ 50 times larger than the signal) forces one 
to restrict the measurements to the "shape-function region", in which the hadronic final 
state has large energy {Ex ~ "2,;,), but only moderate invariant mass (~ rrib^^QCD)- It 
is therefore of great interest to consider differential decay distributions, from which it is 
possible to derive predictions for partial decay rates with arbitrary cuts. In this context, 
a first important set of results was obtained in [13], where it is possible to find analytic 
expressions for the NLO triple-differential distribution of the semileptonic B I v 

decay together with several double and single differential distributions for the same process. 
The resummation of threshold logarithms to next-to-leading approximation in the b ^ u 
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transition was considered in [14]. Higher order contributions to B ^ X^ly decays were 
considered in [15-18] and, very recently, the full NNLO QCD corrections to the partonic 
process b ^ civ were obtained in [19]. Since the OPE applies only to sufficiently inclusive 
quantities, different frameworks were developed in order to account for effects due to cuts 
on the kinematic space [20-26]. In particular, in the shape-function region. Soft Collinear 
Effective Theory (SCET) provides an appropriate framework for the evaluation of the 
triple-differential distribution of the inclusive scmilcptonic decay B Xu I v. The NLO 
analysis of the latter process within the SCET approach is presented in [22,23]. At NNLO, 
the situation is more complicated, but the jet and soft functions are known to 0{ag) in 
perturbation theory [27,28]. The only missing piece is the hard function, which can be 
obtained from the two-loop QCD corrections to the decay of a 6-quark into a u-quark and 
an off-shell W^-boson [29]. On the other hand, these virtual corrections can be considered 
as a first step towards an exact evaluation of the NNLO QCD corrections to the heavy-to- 
light quark transition. To complete the latter calculation, it is also necessary to take into 
account the real emission. 

In this work we focus on the calculation of the two-loop QCD corrections to the decay 
process h u W* . We provide an analytic expression for the three independent vertex form 
factors characterizing the coupling of the quark current with the charged weak boson. These 
form factors arc evaluated by employing a set of techniques which are by now standard in 
multiloop calculations (see for instance [30]). We generate the relevant Feynman diagrams 
with QGRAF [31]. The form factors are extracted directly from the Feynman diagrams 
by means of projector operators. The whole calculation is carried out in Dimensional 
Regularization (DR); UV and IR (soft and collinear) divergencies appear as poles in {d—4), 
where d is the space-time dimension. Since we work in DR, a prescription for handling 
the matrix 75 in d-dimensions must be chosen. We employed a 75 which anticommutes 
with 7^ in c?-dimensions. This prescription is appropriate for the case under study, since 
it is known that the diagrams that we consider fulfill a canonical (non-anomalous) Ward 
identity. After applying the projectors, the contribution of individual Feynman diagrams 
to the form factors is given by a combination of dimcnsionally regularized scalar integrals. 
These integrals are related to a small set of master integrals (Mis) by means of the Laporta 
algorithm [32] . The Mis are evaluated by employing the Differential Equations method [33] 
and they are expressed as Laurent series in {d — 4). The coefficients of the series are 
given in terms of Harmonic Polylogarithms (HPLs) [34] of a single dimensionless variable 
y = q^/ml = —Mi/ml, where is the squared momentum carried by the VF-boson and 
Ml is the lepton pair invariant mass. Since y is negative in the physical region (— 1 < y < 0) 
we perform an analytic continuation y —x — iO~^ , where now x = Mf/ml, < X < 1. 
The analytic continuation is indeed completely trivial, since all the HPLs appearing in the 
result are real for — 1 < y < 0. The form factors found with the above procedure still 
contain UV and IR divergencies. It is possible to get rid of the UV divergencies by means 
of the renormalization procedure. We renormalize the form factors in a mixed scheme: 
the heavy- and light-quark wave functions and the heavy-quark mass arc renormalized 
in the on-shell (OS) scheme, while the strong coupling constant is renormalized in MS 
scheme. The results shown in this paper contain IR divergencies. In order to cancel them. 
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it is necessary to combine these results with the appropriate jet and soft functions [22]. 
Analogously, one can add the exact real emission and consider physical observables which 
are sufficiently inclusive with respect to the hard and soft radiation. 

The paper is structured as follows. In section |2|, we introduce the Feynman diagrams 
involved in the calculation and we discuss their structure in terms of form factors. In 
section |3|, it is possible to find the details of the UV renormalization procedure. In sections 
Q and H, we collect the analytic expressions of the UV renormalized one- and two-loop 
QCD corrections to the form factors, respectively. The expressions of the bare form factors 
as well as the contributions of the individual diagrams to the form factors can be found 
in [35]. In section |6|, we discuss the Ward identity relevant for the b — > uW* decay and we 
prove that our form factors satisfy it. We also provide the analytic expression of the one- 
and two-loop QCD corrections to the scalar vertex in which the pseudo-Goldstone boson 
couples to the quarks, since it enters in the Ward identity fulfilled by the bWu vertex. Our 
conclusions can be found in section |^. Finally, in appendix |A| we collect the set of Mis 
employed in the calculation. 

2. Feynman Diagrams and Form Factors 

We consider the decay process b — > uW* ulvi. The bottom quark of mass carries a 
momentum P and decays in an up-quark (considered as massless) which carries momentum 
p and a VF-boson of momentum q = P—p. Subsequently, the VF-boson decays in the pair Wi 
of squared invariant mass = -q^. The mass-shell conditions are such that P = —mr 
and = 0. The Feynman diagrams contributing to the two- loop QCD corrections to the 
decay process b ^ u W* are shown in Fig. |l|. The most general vertex correction in the 
Standard Model can be described in terms of six form factors Fi and Gi {i = 1, 2, 3): 

y^(P,p) = F,{q')j'^ + -l-F2(g2)^M-g^ + + Giiq^)r75 + ^G^iq^Hq'' 

+ 7^G;iq'h,r , (2.1) 



where = P^ +Pfj,- The spinors u{p) and u{P), multiplying Eq. (|2.1[) from left and right, 
respectively, are not written down explicitly. We define a^'^ = —i/2[y^,^'^]. Since the 
n-quark is taken as massless, only three of the above form factors are independent. By 
replacing 

u{p) h^nl u{P)qu = 2imbu{ph''u{P) - 2qi'u{p)u{P) , (2.2) 
in Eq. (|2.1|), we find the following relations among the form factors in Eq. (|2.1|): 



F2 — — G3 , F3 — —G2 , Fi + -F2 — Gi . (2.3) 

Consequently, using the definitions^ = (1 + 75)72 and Pr = (1 - 75)72, we can rewrite 
the vertex structure as follows: 

V^{P,p) = 2Gi{q^)rPL - —G2{q^)PRq^ - —G3{q^)PB.r ■ (2-4) 

mb nib 



^We employ the notation and conventions of [44]. In particular, in our notation 75 = —75"^, where 75^ 
is the matrix commonly employed in the Bjorken-Drell notation. 
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The form factors are expanded in powers of a^: 



Gi 



2^/2 



ub 



V vr 



(2/) 

The purpose of the present work is to evaluate G- where i = 1,2,3. Vub represents the 
CKM matrix element and the weak interaction coupling constant. With the normaliza- 
tion chosen in Eq. ( |2.5| ), one finds 



Gf^ + O 



a'; 



vr- 



(2.5) 



G'; 



and G 



G 



0. 



(2.6) 



The contribution of the virtual two-loop corrections to the hadronic tensor of [13] can be 
obtained by means of the following translation rules: 



(2/,vir.) 



2 2Gf) + (Gl- 
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(2/,vir.) 



(2/,vir.) 



1 



mt 



G','+ G 



(20, ^(2/), 



G\"'[G)^''+ G 



1 — x 



G. 



{II) 



GfA- 



(2.7) 
(2.8) 
,(2.9) 
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.(10^ 



(2.10) 



Consequently, the hard functions Hij, as defined in [22], can be extracted using the relation 
between Eqs. (16) and (17) of the same article. Note that the hadronic form factor W2 
does not receive contributions from two-loop virtual corrections. 

We write our analytic results in terms of the dimensionless variable 

x = --^ = ^ , 0<ar<l. (2.11) 

In writing or results, we employ the Harmonic Polylogarithms as defined in [34]; on top 
of the canonical weights, we introduce the weights —2 and 2, arising from the integrating 
factors l/(x+2) and 1/(2 — x), respectively. Actually, only two HPLs containing the weight 
2 appear in the final result; they are 

H{2; x) = - In (1 - I) = -i/(-l, 1 - x) + ln(2) , 



H{2, 1, 1; x) = ^^^(^ (1 - 



= ln(l - xf ln(2 - x) - ln(l - x)Li2(-l + x) + Li3(-1 + x) + (3) , 

= -iJ(-l, 0,0,1 -x) + ^C(3). (2.12) 

For convenience, all the results of the paper, including the renormalized and bare 
form factors, as well as the contributions of individuals diagrams, are collected in the file 
SemilepFF . txt [35] included in the arXiv submission of the present work. 

3. UV Renormalization 

The UV renormalization is performed by subtracting the one-loop sub-divergencies and the 
two-loop over-all divergencies. We renormalize the heavy- and light-quark wave functions 
and heavy-quark mass in the on-shell (OS) scheme, while the coupling constant as is 
renormalized in the MS scheme. 

Neglecting for the time being mass renormalization, the bare and renormalized form 
factors satisfy the relation 

G = 4«46Gbare(«r ) > (3-1) 
where in the functions G wc dropped the subscript i = 1,2, 3. 

The perturbative expansion of the various quantities in the equation above is 



G = ^Vub (g(00 + aG(i') + a^G^^^) + O 1 

Z2,u = 1 + aoSzi]l^ + al5zf^ + O (ag) , 
^2,6 = 1 + ao^Z^y + al5Zf^ + O (ag) , 
ao = a (1 + a(54"^ + a^J^^') + O (a^)) , (3.2) 
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where we defined 



ao 



a. 



a 



(3.3) 



(3.4) 



vr vr 
Therefore, the one-loop renormahzed amphtude is given by 

~ '^bare + 2 ' 

where we aheady took into account the fact that 6Z^^ = in the on-shell scheme. The 
two-loop renormahzed amplitude reads instead 

2 ^'" 2 2 

.(10 



G(20 = ) + 
bare 



G 



bare ' 



(3.5) 



To account for mass renormalization, it is sufficient to add the contribution of the 
counter term diagram in Fig. ^ to the r. h. s. of the equation above. 
The renormalization constants are the following: 

1 



6z'^%{d) 



brriQ^ (d, m, 



-C{d) 



d-A 



^lcA + \TR{Ni + Nh) 
3 



mC{d) 



sz^^^dA 



C{d) 



m? J 

2 \ (4-rf)/2 



2 (d-4)((i-3) 

Cf jd - 1) 
~Y{d-4){d- 3) ' 



2 \ 4.-d 



1 



4:~d 



c-\d) ( ^ ) % 



2(d - 4) 24 
Cpfi + Ca/2 + + \NhfA 



where jj. is the renormalization scale and the constants /i , • • • , /4 are [36] 



/i 

/2 

/3 
/4 



9 



51 



8(d - 4)2 32((i - 4) 

11 101 

+ 



128 ~ 2*^^ ~ ^ -* 



13 
16 



^2 + 0(d-4) 



(3.6) 
(3.7) 
(3.8) 
(3.9) 
,(3.10) 

(3.11) 



8(d- 4)2 32(d-4) 128 4 
1 9 59 7r2 

2{d - 4)2 " 8(d - 4) + 32 ^ 12 ^ ^ ^ ^ ' 



+ ?C3 - y ln(2) + + o(ci - 4) , (3.12) 



1 



19 



1139 7r2 

{d - 4)2 24(d - 4) ^ W - y + '^(^ - 4) • 



The factor C{d) is 



G(d) = (47r)(^-'^)/2pJ^, 



(3.13) 

(3.14) 
(3.15) 

(3.16) 



After UV renormalization, the vertex form factors still contain poles in \/{d — A), which 
are associated to soft and collinear singularities. 
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Figure 2: Mass-renormalization counter-term. 



4. One-Loop Form Factors 

In this section we collect the analytic expression of the one-loop renormalized form factors 
defined in Eq. (2.5). In the formulas below, Cp = (iV^ — l)/2Nc is the Casimir operator of 
the fundamental representation of SU{Nc), where Nc is the number of colors (in the SM 
iVc = 3). 

The form factor G^^^^ is given by 

2 1 



GY'^ = G{d) GfY. G^l'^'\d -Ay + O {{d - 4)2) , (4.1) 



m , 

i=-2 



where the first four coefficients in the expansion in {d — 4) are 
Gi^''"^) = -1, 

Gr^^ = ^ + ^(i;x), 

G^l'''^ =-l+ x) - \h{^, 1; x) - H{1, l;x), 



Cr^' = 1- ^-^H{1- x) - 1; x) - ^-^H{1, 1; x) + \h{0, 0, 1; x) 

+^H{0,l,l;x) + ^Hil,0,l;x) + H{l,l,l;x) . (4.2) 
The form factor g'^''^ is 



G^"^ = C(d) ( £^ ) Cj.J2G^;^''\d-4y + 0{{d-4f) , (4.3) 

with 



i=0 



Finally, the form factor Gg^''' is 

Gf =G(d)(4) ' G^^G(^''^)(d-4)^ + 0((d-4)2) , (4.5) 
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where 

G'i'''^ = -Hil; x) + ^H{1, 1; x) + -^H{Q, 1; x) . (4.6) 
X Ix 4x 

Note that the IR poles of the form factor G^^'^ exponentiate. This means that from 
the 0(c?g) expansion of the form factor 

^ = exp{Gi^')}, (4.7) 

we can predict exactly the l/{d — 4)^ and \/{d — 4)'^ poles of the part of the two- loop 
form factor G\ (Eqs. ( |5.3|j5.4| ) below). Moreover, exponentiating also the finite part of 
G^^''' , the double pole of Eq. (|5.5| ) is exactly recovered. 



5. Two-Loop Form Factors 



In this section we collect the analytic expression of the two-loop renormalized form factors 
defined in Eq. (|2.5D . In the expressions below, Ca is the Casimir operator of the adjoint 
representation of SU{Nc), Ca = ^c, Tji is the normalization factor of the color matrices, 
Tji = 1/2, A'^; is the number of massless quarks in the theory, and A'^^ is the number of 
quarks of mass mi,. Therefore, for the decay b uW* in the SM, A'^; = 4, and A'^^ = 1. 
In the finite part of the form factors given below, the constant /C is a rational number. Its 
numerical value is /C = 3.32812ib0.00002, and its analytical value is likely to be /C = 213/64. 
We observe that the following formulas involve HPLs of argument x and of maximal weight 
4. If desired, the HPLs appearing in the equations below can all be rewritten in terms of 
product of Nielsen Polylogarithms of more complicated argument. Because of the chosen 
renormalization scheme, our results depend on the renormalization scale In the formulas 
below, we employ the following notation: 



The form factor G^^''* can be written as 



Gf ^ = C\d) 



4-d 



Cf Gf''\d - 4)* + O (d - 4) , 



(5.1) 



(5.2) 



where the coefficient of the expansion in {d — 4) (up to the finite term) are 

(5.3) 

(5.4) 
49-66L/,+9C(2) 



^ 2 
Gf = Cf 



---H{l-x) 



Ca^^ + TrNi\ 



G 



(21,-2) 
1 



Cf 



7"? 1 _Sr 1 



+ Ca 



72 
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+ TrNi 









3 



(5.5) 



3(71 + 8C(2) - 16C(3)) ^ 13^^^,^ ^ i^^^^^ 



64 



16x 



8x 



+ — H{l,l;x)--H{0,0,l;x)--H{0,l,l;x)-H{l,0,l;x)-4H{l,l,l;x) 

4x 4 2 



1549+1980L^-396L2+972C(2)-1188C(3) ^ 67+66L^-18C(2) 



1728 



72 



-H{l;x) 



-125 - 180L/, + 36L2 - 108C(2) 5 + 6L^ 



432 



18 



if(l;x) 



-5L^ + L2-C(2) 



12 



1 



;(6635(x - 1)^ + (80(x - 1)(-59 - 62x + 3^2) - 480 In (2) (10 



1280(x - 1)3 

22x + 5a;2 + 4a;3))C(2) - 16(261 - 64/C + 41a; + 192/Cx + 405x^ - 192/Ca;2 
99x^ + 64/Ca;3)C^(2) + 40(52 - 72a; + Ibx'^ + 14a;3)((3)) + - — ^T2-(C(2) 

(-49 + 251x - 355x2 



-3xC(2) + x2c(2) + 5x-\(2))F(-l; x) + 



1 



32(x - l)2x 



+153x^ + 12C(2) - 160xC(2) - 148x^^(2) + 24x^C(2) - 16xC(3) + 32x2C(3) 
-16x3C(3))iI(l;x) + 3(-2C(2)-2xC(2)+x2C(2))^^,^^^ 



8(x - 1)3 

C(2)-5xC(2)+3x\(2)-x3C(2) 
-I — — H{0, -1; X) + 



(15 - 106x 



2(x-l)3 ^ ' 32(x-l)3x 

+248x2 - 138x^ - igx'' + 72xC(2) + 152x2C(2) + 48x^C(2))-H"(0, 1; x) 



25 - 134x + 59x' 



H{l,l;x) + 



16(x — l)x 
1 + 27x - 9x2 _ 25^3 + I2x^ 

16(x — l)3x 
l-x+21x2-7x3 



1 - 3x + x"' + 5x'' 
2(x - l)2x 



i^(-l,0,l;x) 



, , 5-59x+83x2-56x3+30x^ ^,^^ ^ 
HiO,0,hx) + _ -^p^ ^(0,l,l;x) 



Hil, 0, 1; x) - ^-^Hil, 1, 1; x) - |t^^i/(2, 1, 1; x) 



4x 



8(x - 1)3 



+ 



4(x - lyx 

l-5x + 3x2-x3^^^ ^ 3 + llx-5x2 + 3x3 

-HiO, -1, 0, 1; x) HiO, 0, 0, 1; x) 



(X - 1)3 

9+9x+ 13x2 -3x3 
4(x - 1)3 



8(x - 1)3 

i7(0,0,l,l;x) + — ^H{0,l,0,l;x)+3Hil,0,l,l;x) 

4(x — Ij-^ 



^i?(0,l,l,l;x)-ilf(l,0,0,l;x)+2i/(l,l,0,l;x)+8i?(l,l,l,l;x) 



+Ca 



29700L2 - 447185- 142560L^ - 3960L3 
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-((11880L^(x-l)' 
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+180(-698 + 1338X - 825x2 + 104x3)C(3)) - ^^^-^^^^^-^^^^ H {-1; x) 
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1d(x — Ij"^ 

_ (l-5x + 3x2--x3)C(2) _ ^ ^ (-33 + 260x + 66L^x 

4(x - 1)3 ^ ' 144(x - l)3a; ^ 

-384x2 - 198L^x2 + 273x3 + 198L^x3 - 1 16x^ - 66L^x^ - 18xC (2) + 468x2 ( (2) 

+72x'^C(2))i^(0,l;x) + —j^ — (-39+235x+132L„x-349x2-132L„x2 

144(x — l)x 

-72xC(2) + 72x2C(2))if(l, 1; x) - ^^:^-t^l±^H{-l, 0, 1; x) 
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(5.6) 
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The form factor 0^2^"^ is 



where 



a 



{21,-2) 



G. 



{21-1) 
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{21,0) 
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{21) 



c- 



2 \ 4-d 



Cf Y1 Gf''\d -Ay + 0{d-A), 
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Cf 
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2x2 

9 7(2 - 5x) 
4x 8x2 
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(-310(x - 1)^ + (60 In (2)x(-38 + 58x - 40x2 + llx^ 



-20(x - 1)(10 - 120x - 79x2 _^ 12x^))C(2) + 16x(125 + 103x)C^(2) 

. . 2 (2-9x-5x2 + 3.T-'5-3x'^)C(2) 

-5x(-30 + llOx - 80x2 _^ 27x^)C(3)) + 

1 



16(X - 1)3x2 

M4x2C(2)-' 
2(2 + x)C(2) 



2(x - 1)'It2 

32 + 195x - 397x2 _^ 337^3 _ ^q^^a _^ 24C{2) - 12xC(2) 



Hi-l;x) 



^A4x^C{2)-76x'%{2)+36x\{2))H{l;x) + 3(30 34x+16x2 3x3)C(2) ^^^.^^ 



iJ(0,-l;x) + 



1 



(X - 1)4 ' ' ' ' 16(X - 1)4x2 

^ - (448x2 + 368x3 )C (2)) 
2 - 9x - 5x2 ^ 3^3 _ 3^4 



4(x - 1)4 
(26 - 69x - 68x2 _ gg^a ^ -^^qq^a 



+3x5 - (448x2 +368x3)C(2))ff(0,l;x) + ^ 18x+37x2 + 172x3 ^^'^V (l,l;x) 



8(x - 1)2x3 



, 30- 34x + 16x2 - 3x3 

(.-1)3.= «(-!.»■ 1;-) + 40— [,3 «P.M;x) 

2 - 27x + 4x2 - 48x3 + 66x4 _ 5 7(2 - 3x) , , , 
8(,^ _ 1)4,^2 ^(0' 0' 1; ^) + ^2^^^^^' ^' ^' ^) 

10 -39x+ 234x2 -276x3 +86x4- 24x5 4(2 + a;)„, , 

H{0, 1, 1; x) + / ,' h{0, -1, 0, 1; x) 



(X - 1)^ 



4(X - 1)4x2 

2- 13x - 22x^ - 12x3 ^^^^ _^ ^i±^i/(0,0,0, l;x 



2(X - 1)3x2 



+ 



6(4 + 3x) 



if(0, 0,1,1; x)- 



3(4 + 3x) 



(x - 1)^ 
^(0,l,0,l;x) 



+C^ 



(X - 1)4 

(1320L^(x - 1)^ + 20(x - l)3(-269 + 242x) 



1440(x - l)4x 

+(540(x - 1)(4 + 50x + 5x2 ^ g^a^ ^ ^^q^^ (2)x(38 - 58x + 40x2 - llx3))C(2) 
+36x(364 + 317x + 108x2)^2(2) + 45x(-30 + llOx - 80x2 ^ 27x3)C(3)) 



(2-9x-5x2+3x3-3x4)C(2) 



H{-l;x) + 



r(406 + 132L^-2067x 



4(x- 1)3x2 ^'"'^ ' 144(x- 1)3x2^""" ' ^ 

- 594L^x + 3603x2+ qqql^x'^_ 2629x3- 726L^x3+ 687x4 + 198L^x4 - ( 144 - gg4^ 

-1224x2-1242x3 + 54x4)C(2))/f(l; x) - 3(30 - 34x^+16x2 - 3x3)C(2) 
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(x-l)4 ^'"^ ■ 24(a;-l)4a;2 

j^- (240x^+2100;^ + : 
2-9x-5x^+3x^-3x'^ 



(-22 + 145a; - SOOx^ + 109a;^ + 59a;^ 



+9a;5-(240x2+210a;3 + 72x^)C(2))ii'(0,l;x) - ^^-J^^p:^^—^Hil,l;x) 



24(x - 1)2^2 

rr.,^, , , 4+ 76a;- 96x2 + 16x^-93;^ ^,^ 
F(-l, 0, 1; x) + 8(x - l)4x ^ 



2(x - 1)3x2 

8 + 58x - 24x2 - 36x3 + 3x^ , , 4 ^_ 24x + 11x2 + Sx^ ^^^^ ^ ^ ^ 

H{0, 1, 1; x) + H{1, 0, 1; x) 



(x — l)3x 



8(x - l)4x 

30 - 34x + 16x2 - 3x3 24+17x + 4x2 ^ 
8(^31)4 H{2, 1, 1; x) + ^(^31)^^(0' 0' 0' ^) 

.^^^±^i/(0, -1, 0, 1; x) + ^±^^±^(2i/(0, 0, 1, 1; x) - i/(0, 1, 0, 1; x)) 



(x - 1)4 ^ ' ' ' ' • ^ ' 4(,T _ 1)^ 

19+6L^ 26-51x+6(2-3x)L^ 



+TrNi 
+TRNh 



18x 
1 



36x2 



2 — 3x 

H{1; x)+-^{H{0,l; x) + 2iJ(l,l; x)) 



18(x - l)4x 



(-19 - 6L^ - 164x + 24L^x + 393x2 _ 3g^^^2 _ 2ig^3 



+24L^x^+8x''-6L^x4+(252x - 300x2+84x^-36x'')C(2) - (72x + 36x2)C(3)) 



+ 



26 - 223x - 124x2 - 51x3 + (12 - 42x + 48x2 - 18x3)L^ 

36(x — 1)2x2 

2 - 9x - 21x2 _ ;l3x3 - 3x4 2(2 + x' 

-i?(0, l;x) + 



6(x — 1)3x2 
.{20 



H{l;x) 



{x-iy 



-H{0,0,l;x) 



(5.10) 



The form factor Gg ' can be written as 



G 



(20 



C\d) 



i-d 



Cf Gf''^ {d -Ay + 0{d-A) , 



with 

Gf = Cf 



2x 



H{l;x) 



G 



(21,-1) _ 



G. 



(21,0) 



Cf 
Cf 



-^H{l;x) - ^{HiO, l;x) + 6iJ(l, l;x)) 



(5.11) 

(5.12) 
(5.13) 



1 



-((-60(x - l)x(57+2x) -60 In (2)(-22+30x-26x2+9x3))C(2) 



80(x - 1)4 

^8(9+49x+18x2)C2(2)-5(22-90x+82x2-41x^)C(3))- 

+161x2 - 55x^ - (12 - 364x - 436x2 - 28x2)C(2))i?(l; x) 
, (l + 5x + 5x2 + x3)C(2) ^^ ^ ^ , 3(-14 + 6x-2x2 + x3)C(2) ^^^_ 

+ -.xQ V -'■)^j+ J/ -|\4 I1[Z,X 



16(x— l)3x 



(51-157X 



2(x - l)3x 



6xC(2) 
(x - 1)4 



H{0,-l;x) + 



1 



16(x— l)4x 



4(x - 1)4 

(11 - 72x+ 366x2 _ 260x-^ - 45x^ + 96xC(2) 



+528x2C(2) + 192x^C(2))i?(0, 1; x) - ^ + 36x + 8x2 24x3 3^4 ^_ 

8(x — Ij^x 
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4(x - l)4x 
14 - 6x + - x^ 



2(x - l)3x 



4(x - 1)4 
2+3x+4x2 



, 2(2 + 15x + 4x2) , 
H{2, 1, 1; x) - ^- ^-HiO, 0, 1, 1; x) 



(x - 1)4 

12x 2+15x+4x2 
4 iJ(0,0,0,l; x)- -^—^ if (0,-1,0,1; x) + ^—^^if (0,1,0,1; x) 



(x-iY 



160(x-l)4 



(60(x-l)^ + (-20(x-l)(56 + 105x + 40x2) + 60 In (2) (-22 



+30x - 26x^ + 9x'^))C(2) - 12(18 + 137x + 108x^)^(2) - 5(-22 + 90x - 82x 



+41x3)C(3))-(^ + ^^ + ^^\-t^'^^(^^/i(-l; ^ ' 



■,x + 



-(335 + 66L^ 



4(x-l)3x ^'""^ ' 144(x - l)3x^""" ' """^ 

-843x - 198L^x + GSlx^ + IDSL^x^ - 173^3 - 66L^x^ - (72 + 468x + 2466x2 
+126x3)C(2))ii(l;x) - 3(-14C(2)+6xC(2)-2x2C(2) + x3C(2))^^^^^^ 



3xC(2) 



HiO,-l;x) - 



(x-l)4 ^ ' ' ^ 24(x-l)4x 
270x2 + 216x3)C(2 
1 + 5.x + 5x2 _^ ^3 



8(x - 1)4 

(11 + 22x - 171x2 + 79x3 _^ 59^4 _ ^gg^ 



+270x2 + 216x3)C(2))i/(0, 1; x) - 20 + 40x^^66x2 3x3 ^^^^ ^_ 



2(x - l)3x 
22 + 56x - 62x2 - 7x3 

8(x - 1)4 
14 - 6x + 2x2 - x3 



TTf ^ n^ ^ 13 + 82x + 58x2 
iJ(-l,0,l;x) — -2 ii(l,l;x) 



iJ(0,l,l;x) 



24(x - l)2x 
10 + 23x + 9x2 
" 8(x - 1)3 



i?(l,0,l;x) 



8(x - 1)4 



Hi2, 1, 1; x) - ^±^l^J^ii(0, 0, 0, 1; x) 
4(x — 1)^ 



6x 2 + 7x + 12x2 
+7 —H{Q, -1, 0, 1; x) ___^(2ii(0, 0, 1, 1; x) - H{{), 1, 0, 1; x)) 



(X - 1)4 

(25 + 6L^) 



4(x - 1)4 



+TrNi 



+TRNh 



ii(l;x) + — (/i(0,l;x)+2i7(l,l;x)) 
(54 - 91x + 20x2 + 17x3 - (36 + 12x - 52x2 + 4x3)C(2) + 36xC(3)) 



25+322x+25x2 + (6- 12x+6x2)L^ 



6x 



36(x — l)2x 



-ii(0,0,l;x) 



6(x - 1)4 

H{l;x) 



l + 21x+21x2+x3 
6(x — l)3x 



■HiO,l;x) 



(5.14) 



We checked our results for the form factors G\ (x) f = 1,2,3 against the calculation 
of Martin Beneke, Tobias Huber, and Xin-Quing Li [37] and we found complete analytical 
agreement. 
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6. Ward Identities 



We explicitly checked that the UV renormalized form factors satisfy the on-shell Ward 
identity^ 





(6.1) 



where (f) is the charged pseudo-Goldstone boson and the gray circles represents the sum of 
all two-loop one-particle-irreducible QCD corrections to the vertices. The Lorentz index 
associated to the VF-boson is saturated by the boson momentum g^. In order to satisfy 
the relation in Eq. (|6.1| ) it is necessary to renormalize also the factor mi, appearing in the 
tree-level (p'^ub coupling. The relevant NNLO mass counter term can be found in [36]. 

The two-loop corrections to the scalar coupling of the pseudo-Goldstone boson to quark 
can be absorbed in a single form-factor S, defined as follows 




' Mw 



S{q^)u{p){l-^,)u{P). 



(6.2) 



The UV renormalized form factor has the following perturbative expansion in a^: 



S 



2^/2 



Vuh 



+ 



TT 



^(10 



+ 



vr 



2 /a? 



with 5(°') = 1. 

The one-loop form factor 5*^^'^ is given by 



(6.3) 



/ ,,2 \ (4-d)/2 1 

SW = C{d)l^] S^^^''\d-4:y + (D{id-4:)^) . (6.4) 

After UV renormalization (including the renormalization of the Yukawa (j)^ub coupling), 
the coefficients of the expansion in (d — 4) are 

-2) 



-1, 



-1 - —H(l;x) 
2x ^ ^ 



H{l,l;x)--H{0,l;x), 



1 + ^^H(l; x) + ^H{1, 1; x) + ^H{0, 1; x) + H{1, 1, 1; x) 
+ \h{1, 0, 1; x) + ^H{0, 1, 1; x) + ^H{0, 0, 1; x). 



(6.5) 
(6.6) 

(6.7) 
(6.8) 



^It can be proved that the Ward identity is fulfilled already at the level of master integrals, irrespectively 
on the analytic expression of the Mis themselves. 
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The two-loop form factor 5^^'^ is given by 



^ 2 \ 



(6.9) 



where the coefficient of the expansion in {d — 4) are 



-H(l-x) 



cj7 2 4- 1 



49 + 9C(2) 



12 ^ 12 ^ ' 



3(47 + 8C(2) - 16C(3)) 7-|-10x 



+ TRNh 



72 
1 



(6.10) 
(6.11) 

(6.12) 



64 



8x 



if(l;x)-^|^iJ(0,l;x) 



6-l-5x 1 3 

■^^H{1, 1; x) - -//(0, 0, 1; x) - -H{0, 1, 1; x) - H{1, 0, 1; x) - 4/7(1, 1, 1; x) 



(1549+1980L^-396L2+972C(2)-1188C(3)) 67+66L,,-18C(2) 



+TrNi 



1728 

-2 



72 



Hil-x) 



- 125 - 180L^ + 36L^ - 108C(2) 5 + 



+TRNh 

831 
256 



432 

-5L^ + L^-C(2) 
12 



18 



Hil-x) 



--L^H(l-x) 



(74 - 96x + 13a;2)C(3) (1 + 2x + x2)C(2) 



+ 



5 ' ' 32(,T - 1)2 2(x - l)x 

17 + 8x - 25x2 - (12 - 78x + 30x2)C(2) - (8x - 8x2)C(3) 

16(x — l)x 

-8x''^ + 12xC(2) + 24x2c(2))i7(0, 1; x 
1 + 2x + x2 



9(5 - 34x + 11x2) 
80(x - 1)2 

H{-l-x) 



(x — l)x 
14 + 5x 



16(x - l)2x 

4-x + 5x2 
^^H{l^■x) 

(-1, 0, 1; x) - 2 -13a: + 20x2 -3x3 ^^^^ ^_ 



H{l-x) 

(11 + 12x - 15x2 



16(x — l)2x 



4x 

2 + 5x - 4x2 



mill ^ ^ 10-^^^+20^! + ^ P7^n 1 1 ^ 
H(l, 1, 1; x) H — — 2 ii (0, 1, 1; x) 



8(x - l)2x 

SfA — Ax + x'^) 

H{l,0,l;x)-\- -^i/(2,l,l;x)-/7(0,-l,0,l;x) 

8(x — 1)"^ 



4(x - l)x 

5 _ 2x + 3x2 ^^^^ 0, 0, 1; x) + i— l^^^±^iJ(0, 0, 1, 1; x) + Ji? (0, 1, 1, 1; x) 



8(x - 1)2 



4(x - 1)2 
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4 — 4x + 3x'^ 1 

+—U — T\^^(0' 1' 0' 1; ^) - o^(i' 0' 0' 1; + 3^^(i' 0' 1' 1; ^) 

4(a; — Ij^ 2 



+2/7(1, 1, 0, 1; x) +8F(1, 1,1,1; x) 

( 

5 y 



+ Ca 



+ 



/IIL^ 31n(2)(x - 4)(7x - 8) 1067 + 49x\ 
~96 16(x - 1)2 576(x - 1) J 



54589 IIL^ 55L2 IIL^ 
20736 12~ ^ T92 288~ 

179 - 250x + 125x2 



C(2) 



2K\ ^ (896 - 1324x + 347x2)C(3) (1 + 2x + x2)C(2) 



576(x - 1)2 



4(x - l)x 



160(x - 1)2 
H{-l;x) 



+ 



1 



;(708 + 198L^ - 466x - 198L^x - 99L2x - 242x2 + 99LIx^ 



432(x - l)x 

-216C(2) + 738xC(2) - 684x^^(2) - 378xC(3) + 378x2C(3))F(l; x) 

+66L^x-211x2-132L^x2 + 89x^+66L^x^-(126x-144x2 + 72x^)C(2))F(0,l;x) 

78 + 223x + 132L^x-72xC(2)^,^ ^ , 1 + 2x + x2 

/i(l,l;x) — — — ii(— l,0,l;x) 



144x 
40 - 56x + 7x2 



11 



48(x-l)2 ^(0,0,l;x)--/7(l,l,l;x) 



2(x - l)x 

44 - 88x + 53x2 



(x - 1)2 



^(0,1, l;x) 



1 



1 



+ -H{l,0,0,l;x) 



TrNi 



4(x - 1)2 
1 



F(0,0,l,l;x) + 



1 



8(x- 1): 



rF(0,l,0,l;x) 



5184 



(3893 + 1728L^ - 540L^ + 72Lj^ + 3420C(2) 



-216L^C(2) + 720C(3)) + 



48 + 18L^ + 28x - 9L2x + 45xC(2) 
108x 



Hil;x) 



+ 



3 + 5x + 3L„x , 3 + 5x + 3LuX , , 1 
'^-F (0, l;x) H 



18x 



^i7(l,l;x) + -/7(0,0,l;x) 



+Ih{0, 1, 1; x) + ^H{1, 0, 1; x) + '^H{1, 1, 1; x) 



+ TnNh 



hi _ ^3: + :S 

3 48 72 



+- 



11407 - 17630X + 8527x2 /L. 409 - 747x + 651x2 - 185x3 



2592(x - 1)^ 



24 



144(x - 1)3 



C(2) 



703) 
36 ' 108(x-l)2x 



+ 



(48 + 18L^. + 104x - 36L^x - 9L^x - 112x2 _^ 18L^.x 



+ 18L2x2 + 56x3 _ g^2^3 ^ gx(^(^2) - 18x2c(2) + 9x3C(2))F(l; x) 

(3+14x+3L^x-9L^x2 - 6x3 _^ g^^^s _^ 5^4 _ sl^x'^-jh{0, 1; x) 



18(x - l)3x 



+^H{l,l;x) + ^H{0,0,l;x) 



(6.13) 



When written in terms of form factors, the Ward identity in Eq. (|6.1| ) reads as follows: 
2Gf^^ (x) + xGf ^ (x) + Gf ^ (x) - 25(2') (x) = . (6.14) 
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It can be checked that the form factors presented in this paper fulfill Eq. ( 3.14 ). 



7. Conclusions 

In this paper, we presented analytic expressions for the two-loop QCD corrections to the 
decay process 6 — > u W* ulv. This process is important for the precise determination 
of the CKM matrix element V^^ and, therefore, for the study of flavor and CP violation 
within and beyond the Standard Model of fundamental interactions. 

The Lorentz structure of the process is parametrized in terms of three form factors, 
whose analytic expression are given in the form of a Laurent series of — 4), where d is 
the space-time dimension. The coefficients of the series are expressed in the well known 
functional basis of harmonic polylogarithms of a single dimensionless variable. The result 
can be used in a SCET framework, after combining it with the jet and soft functions 
already known in the literature, for a phenomenological determination of \ Vub\- The results 
presented here are the first step towards a complete determination of the NNLO QCD 
corrections to the heavy-to-light quark transition. 
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A. The Master Integrals 

In this Appendix we collect the analytic expressions of the Master Integrals for the Feynman 
diagrams of Fig. ^. We provide only eight of them, since the other Mis can be found 
in [41,42]. It must be pointed out that the Mis (a)-(f) in Fig. ^ were already calculated 
in [43]. We checked the analytic expressions that we obtained against the results in [43] 
and we found complete agreement. Moreover, all the Mis were checked by comparing their 
numerical value to the results obtained by direct numerical integration with the sector 
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(a) 



(b) 



(c) 



(d) 







Figure 3: Master Integrals needed for the Two-loop QCD corrections. Tiiick lines represent 
massive particles, thin lines represent massless ones. 



decomposition method. The numerical integration was carried out by using the package 
FIESTA (see [40]). The checks were done for several values of the variable y. 

The explicit expression of the Mis depends on the chosen normalization of the inte- 
gration measure. The integration on the loop momenta is normalized as follows 



1 



C{d) 



(rf-4) 
2 



(4^2 



(d-2) ' 



(A.l) 



where C{d) is defined in Eq. ( |3.16| ). In Eq. ( [A.l ) fi stands for the 't Hooft mass of dimen- 
sional regularization. The integration measure in Eq. ( A.l[ ) is chosen in such a way that 
the one-loop massive tadpole becomes 

^ , = , . ^ . (A.2) 



A;2 + m2 {d-2){d-4)' 

In the expressions below, /C is a rational number (its numerical value is /C = 3.32812 it 
0.00002 ~ 213/64), while 04 = Li4(l/2) = 0.51747906.... C(2) and C(3) are the Riemann C 
function evaluated in 2 and 3 respectively: C(2) = 1.6449341..., C(3) = 1.2020569.... 

The expressions of the Mis are the following. 




m^(l + y)2 



Y^Aiid- Ay + 0{d - 4) , 



(A.3) 



(A.4) 



= (A.5) 
b 

= -^C(2) + ^i/(-l,-l;y), (A.6) 

8Q 7 2 

= _C(3)--C(2M-l;y) + -i/(-l,-l,-l;y), (A.7) 

= -^C'(2)/C + ^C(3)/?(-l; y) - ^Q{2)H{-1, -l-y) + \h{-1, -1, -1, -1; y) 

_li7(_l,0,0,-l;y). (A.8) 



= . i?o + 0(rf-4), (A.9) 

^0 = ^^'(2) + 4c(2)^/'(0, -1; y) - ^H{Q, -1, 0, -1; y) + \h{{), 0, -1, -l-y) 
loU lb lb 8 

--^F(0,0,0,-l;y). (A.IO) 



C-i = ^C(2) , (A.12) 

Co = 7?7C(3) + 4c(2)^(-i; y) + 4^(0' -1' -1; y) - 4^(0' 0' -1; y) > (a-i3) 



= "^"^'^^^ + ^C(3M-i;y) - ^C(2M0,-i;y) + ^C(2M-i, -i;y) 

+ ii/(-l, 0, -1, -1; y) - \h{-1, 0, 0, -l-y) + -^//(O, -1, -1, -1; y) 
8 8 lb 

-4^^(0, -1, 0, -1; y) - |-F(0, 0, -1, -1; y) + |-//(0, 0, 0, -1; y) . (A.14) 



(A.15) 



i=-2 



■2 — 77 



i?-i = -- + ^^(-i;y). 

lb 8y ^ ' 



(A.16) 
(A.17) 
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i + - + ^^(-'- -'^ 



+- 



+ 



1 



+ 



64 



C(3) 



^gC(2) ln(2) - Ac(2)^(-2; y) - ^H{-2, -1, -1; y) - ^H{0, -1, -1; y) 

+^H{0,0,-l;y) , (A.18) 
lb 

fi^ 7 S 1 Q 

-64 - 128^^^) - 32^(^) - 32^(^) ^^^^^ + 32^(-^' + 32^(2)^^-^^ ^) 
-^(^(2M-2;y) + A//(o,-l;y) _ ^H{-l,-l;y) - 1f(-2, -1, -1; y) 

-1, -1; y) - Ih{-1, 0, -1; y) - ^H{0, -1, -1; y) + i,/7(0, 0, -1; y) 
16 8 32 16 



+ - 



1 Q 3 1 ^ 

-i/(-l;y) + _C(2M-l;y) - -H{-l,-l;y) + -H{0,-l;y) 



~H{-1, -1, -1; y)- ^H{-1,0, -1; y)- ^//(q, _l, -l-y)+ Ih{0, 0, -1; y) 
16 16 16 8 



+ 



1 



(1 + y) 



1 

96 



QQ 1 

ln^(2) - --C(3) + :^C(2)ln(2) + -^C(2)ln2(2) + — C'(2) - -a. 



5 

128' 



3 

32' 



1 

16" 



640 



^C(3) - ^C(2)ln(2))//(-l;y) + (^C(3) + ^C(2)) ^(-2; y) 

-^C{2)H{0, -1; y) + ^(:(2)H{-1, -2; y) - lc(2)^(-2, -1; y) 

+^i7(-2,-l,-l;y) + 1^(0, -1, -1; y) - 1^(0,0, -l;y) 

-^H{-2, -1, -1, -1; y) + 3^if(-2, -1, 0, -1; y) + ^^^(-1, -2, -1, -1; y) 

+-i:/7(-l,0,-l,-l;y) - -^^^(-1, 0, 0, -1; y) - -|f(0, -1, -1, -1; y) 
16 16 16 



+^H{0, -1, 0, -1; y) + ^H{0, 0, -1, -1; y) - ^F(0, 0, 0, -1; y) 



(A.19) 



1 2 

— V ii;^ (d - 4)* + Oid - 4)3 , 



(A.20) 



(A.21) 



1 



16(2 + y) 



3C(2)+/f(-l,-l;y) 



(A.22) 



-^Hi-1, -1, -1; y) + ^H{-1, 0, -1; y) - ^HiO, -1; y) + ^F(0, -1, -1; y) 



1 



^(0,0, -l;y) 



+ 



1 



+ 



■^C(3) + ^C(2) ln(2) + Ac(2)i?(-2;2/) 



+^iI(-2, -1, + ^H{0, -1, -1; y) - ^H{0, 0, -1; y) 



+ 



1 



(2 + y) 



+-H{-i,o,-i;y) 



(A.23) 



-£^2 — — 

y 



I + T^^(^) + 1^^^) + |c(2)ln(2))i.(-l;y) + (i + |c(2)) //(O, -1; y) 
A + ^C(2))i?(-l,-l;y) + Ac(2)iJ(-l,-2;y) - 0, -1; y) 

+^if(0,0,-l;y) - ^ii'(0,-l,-l;y) + -1, -1; y) 

+ -^iJ(-l,-2,-l,-l;y) - AiJ(-l,-l,-l,-l;y) + l/f(_l,_l,0,-l;y) 
+^i/(-l,0,-l,-l;y) 0,0, -l;y) + ^ii'(0, -1,-1, -l;y) 

-^iJ(0, -1, 0, -1; y) - ^H{0, 0, -1, -1; y) + ^H{0, 0, 0, -1; y) 



+ 



1 



+ 

5 

64' 



— ln^(2) + — C(3) - — C(2) ln(2) - — C(2) ln'^(2) - — C^(2) + -04 
96 ^'' 64^^^ 16^^ 16^^ ^' 640^ ^ ^ 4 



+ (^^^^^ " ^C(2)ln(2))if(-l;y) - (^((3) + ^C(2))i?(-2; y) 
+^C{2)H{0, -1; y) + ^C{2)H{-2, -1; y) - ^(;{2)H{-1, -2; y) 
-^il(0,-l,-l;y) + ^il(0,0,-l;y) - ^i?(-2, -1, -1; y) 
+^iJ(-2,-l,-l,-l;y) - liJ(-2,-l,0,-l;y) - 1f(-1, -2, -1, -1; y) 
-^i/(-l, 0, -1, -1; y) + ^ii-l-l, 0, 0, -1; y) + ^ii-(0, -1, -1, -1; y) 

—HiO, -1, 0,-1; y) - ^H{0, 0, -1, -1; y) + ^H{0, 0, 0, -1; y) 



- 21 - 



+ 



(2 + y) 



7 S 4'i /I ^ 



+^C(2) ln(2)) H{-1- y) - (3^ + ^C(2)) ^(-1, -1; 2/) + ^C(2)i/(-l, -2; y) 

+ _i, _i; y) _ -^^^(-l, 0, -1; y) + ^H{-1, -2, -1, -1; y) 

lb lb lb 

--^/?(-l, -1, -1, -1; y) + _i, 0, -1; y) + |-F(-1, 0, -1, -1; y) 

lb lb 62 



(A.24) 



m^(l + y) 



Y^Fiid- iy + 0{d - 4)2 , 



(A.25) 



F_i = ic(2) + ^^(0,-l;y), 



(A.26) 



^0 = -^^(3) - ^C(2) ln(2) + ^C(2)^(-l; y) - Y^mH{-2; y) - YqH{-2, -1, -1; y) 

+^ii'(-l,0,-l;y) + ^//(0,-l,-l;y) -^F(0,0,-l;y), (A.27) 

1 1 997 1 7 

Fi = -- ln4(2) + -C(2) ln2(2) + -^C'(2) - 704 - ^C(3)//(-l; y) + :i-C(3)F(-2; y) 



96 



16 



640 



16 



32 



+ ^[3F(0,-l;y) +4F(-l,-l;y) -5/f(-2,-l;y)] - F(-2, -1, -1, -1; y) 
61 16 

+-^if(-2, -1,0, -l;y) + \h{-1, -1,0, -l;y) + -H{-1,{), -1, -l;y) 
16 8 4 

-^/f(-l, 0, 0, -1; y) + ^F(0, -1, -1, -1; y) - ^H{Q, -1,0, -1; y) 



--|/f(0, 0, -1, -1; y) + ii7(0, 0, 0, -1; y) 
16 8 



(A.28) 




Y,Gi{d- 4)^ + 0{d - 4)2 , 



i=-2 



(A.29) 



G_2 

Go 
Gi 



-, (A.30) 

-— , (A.31) 

16 ^ ^ 

19 1 1 1 

6^ S ^ 1 ^ 1 



- 22 - 



ii^(0,-l,-l;,) + |i7(0,0,-l;,) + ^ 



+^C(3M-1; y) - 3^C(2MO, -1; y) + ^^^(0, 0, -1; y) + ^//(-l, 0, 0, -1; y) 



~H{0, -1, 0, -1; y) - ]h{0, 0, -1, -1; y) 
lb 4 



(A.33) 




m2(l + y) ; 



^ J,(d-4r + 0(d-4)2 



j=-i 



(A.34) 



J_i = ic(2) + ^^r(0,-l;y), (A.35) 

^0 = -^C(3) + ^C{2)Hi-l;y) + ^i/(-l,0, -l;y) + i//(0, -1, -1; y) 

-^/7(0,0,-l;y), (A.36) 
lb 

23 1 1 1 

+^ii-(-l,0, -1, -l;y) - ^H{-1,0,0, -l;y) + ^i/(0, -1, -1, -l;y) 
4 lb 2 

- Jli-(0, -1, 0, -1; y) - ^H{0, 0, -1, -1; y) + ^i/(0, 0, 0, -1; y) . (A.37) 
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